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We report calculations designed to assess the ultimate precision of an atomic clock based on the 
578 nm 6 ^5*0 — » G^Pq transition in Yb atoms confined in an optical lattice trap. We find that this 
transition has a natural linewidth less than 10 mHz in the odd Yb isotopes, caused by hyperfine 
coupling. The shift in this transition due to the trapping light acting through the lowest order AC 
polarizability is found to become zero at the magic trap wavelength of about 752 nm. The effects 
of Rayleigh scattering, higher-order polarizabilities, vector polarizability, and hyperfine induced 
electronic magnetic moments can all be held below a mHz (about a part in 10^*), except in the case 
of the hyperpolarizability larger shifts due to nearly resonant terms cannot be ruled out without an 
accurate measurement of the magic wavelength. 



PACS numbers: 06.30.Ft, 32.80. Pj, 31.15.Ar 

Optical atomic clocks offer new opportunities for cre- 
ating improved time standards as well as looking for 
changes in fundamental constants over time, measuring 
gravitational red shifts, and timing pulsars. Compared 
with microwave atomic clocks, optical clocks have the 
intrinsic advantage that optical transitions have a much 
higher frequency and potentially much higher line-Q than 
microwave transitions. Moreover, optical frequency comb 
techniques P] now permit different optical frequencies to 
be compared with each other at the 10~^^ level or bet- 
ter, and to be linked to microwave clocks as well. Optical 
transitions in alkaline earth atoms offer remarkable pos- 
sibilities for clocks. In addition to the relatively sharp 
^5*0 intercombination line already in use |3| , there 

occurs the much sharper ^Sq — > ^Pq line in odd isotopes. 
This one-photon transition is forbidden in even isotopes, 
but in odd isotopes acquires a weak El amplitude in- 
duced by the internal hyperfine coupling of the nuclear 
spin. Doppler and recoil shifts can be virtually eliminated 
by confining very cold atoms in an optical lattice trap. 
This lattice will be Stark-free if it is produced by laser 
beams tuned to the magic frequency where the ground 
and excited states undergo the same light shift, leaving 
the clock transition unshifted and relatively insensitive 
to the laser polarization. 

Katori 3] has pointed out these advantages for ^^Sr, 
and the 5 ^5*0 5 "^Pq transition in this isotope has re- 
cently been observed 4J. Also, Sr atoms have been held 
in an optical lattice at the laser wavelength appropri- 
ate for a Stark-free 5 ^5*0 — > 5 ^Pf transition, and this 
transition has been observed free of Doppler and recoil 
shifts . The magic frequency for the Sr clock has been 
evaluated recently in Ref. and it has been measured 



in Ref. In Ref. |^ an estimate of systematic uncer- 
tainties for Sr also has been carried out. 

Ytterbium has two stable odd isotopes, -'^^^Yb and 
i73yi-,^ which also appear to be excellent candidates for 
an atomic standard, using the 6 ^5*0 — > 6 ^Pq transition at 
the convenient wavelength 578 nm. The atoms are read- 
ily trapped into a MOT operating on either the strong 
6 ^5*0 — > 6 ^Pf line or the 6 ^5*0 6 ^P^ intercombina- 
tion line and in the latter case have been cooled to iiK 
temperatures by Sisyphus cooling Also, these iso- 
topes have been successfully confined in an optical dipole 
trap 1^ . In this paper we present a calculation of the nat- 
ural linewidth of the clock transition, which turns out to 
be less than 10 mHz, and also a calculation of the Stark- 
free wavelength of an optical lattice trap for the clock 
transition. This wavelength turns out to be about 752 
nm, reachable with adequate power by a tunable Ti:Sa 
laser. We also estimate the size of the polarizability due 
to higher magnetic dipole and electric quadrupole optical 
moments. In addition, we estimate the Rayleigh and Ra- 
man scattering rates in the optical lattice which can limit 
the coherence lifetime of the clock transition. Finally, we 
compute the small but important hyperfine-induced Zee- 
man shift in the excited state and the vector light shift 
which can cause a small Stark-shift dependence on the 
polarization of the trapping light. 

Our results indicate that a light intensity of 10 
kW/cm^ would create a convenient trap depth of 15 /iK 
at the magic wavelength, while perturbations to the clock 
frequency could be held below the mHz level (10~^* rel- 
ative shift) - with one possible exception. Larger shifts 
due to accidental near resonances in the hyperpolarizabil- 
ity cannot be ruled out without an accurate measurement 
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of the magic wavelength. 

AU the calculations reported in this paper have been 
carried out using the relativistic many-body code de- 
scribed in Refs. [iflliHl3- The employed formalism is 
a combination of configuration-interaction method in the 
valence space with many-body perturbation theory for 
core-polarization effects. The effective core-polarization 
(self-energy) operator is adjusted so that the experi- 
mental energy levels are well reproduced. In addition, 
the dressing of the external electromagnetic field (so- 
called core shielding) is included in the framework of the 
random-phase approximation. In the following we refer 
to this many-body method as CI-hMBPT. For Yb, the 
CI-I-MBPT method has an accuracy of a few per cent 
for electric dipole matrix elements and magnetic-dipole 
hyperfine constants Unless specified otherwise, 

we use atomic units (|e| = h = nif, = f ) throughout the 
paper. 

In the proposed design, the Yb atoms are confined to 
sites of an optical lattice (formed by a standing- wave laser 
field of frequency to and amplitude of electric field £q). To 
the leading order in intensity and the fine-structure con- 
stant, the laser field shifts the clock transition frequency 
Wo by 

5u,{u^) = [allp. {u) - ae%, (c.)} {Eo/2f , (1) 

where (w) is an a.c. electric-dipole polarizability of 
level X 
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FIG. 1: Electric dipole a.c. polarizabilities for 6 (solid 
line) and 6 ^Po (dashed line) states of Yb. The polarizabilities 
are shown as a function of laser frequency ui. 
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with Q^'^'^) being relevant multipole operators fl^. Typ- 
ically, the EI polarizability (|2Jl overwhelms this sum. 
Compared to the EI contribution, the higher-order mul- 
tipole polarizabilities are suppressed by a factor of 
(aw)^'-'"^ for EJ and by a factor of Q;^(aw)^"'"^ for MJ 
multipoles. We verified that at the magic frequency there 
are no resonant contributions for the next-order E2 and 
Ml polarizabihties and we expect a(^2,Mi) < io-6q,(si) 
similar to the case of Sr 0. At the same time we no- 
tice that a core-excited state 4/-^^(^Fi',„)5d5/26s^ J = 5 



We have carried out the calculations of the EI a.c. polar- 
izability using the CI-I-MBPT method. We summed over 
the intermediate states in Eq. ^ using the Dalgarno- 
Lewis-Sternheimer method [13 • The results of the calcu- 
lations for both 6 ^5*0 and 6 ^Pq states are shown in Fig.^ 
The two dynamic polarizabilities intersect at uj* = 0.0606 
a.U. At this "magic" frequency the lowest-order differen- 
tial light-shift, Eq. vanishes. It is worth noting that 
at UJ* the sum (0) for the ground state is dominated by 
the 6s6p state and for the 6 ^Pq level by the 6s7s ^Si 
state. 

In general, for a linear laser polarization the second- 
order light shift of level X can be represented as a sum 
over 2'^-pole polarizabilities a^x^"^ (-^ distinguishes be- 
tween electric, A = 1, and magnetic, A = 0, multipoles) 



SEx = -^j:a^ 
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(C) 
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When the total angular momentum of the level X is equal 
to zero these a.c. polarizabilities are expressed as 



J + 1 2J+1 
J [(2J+1)!!] 
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may become resonant with an excitation from the 6 '^Pq 
level. The relevant M5 polarizability is highly sup- 
pressed, and we anticipate that the magic frequency will 
be only slightly shifted by the presence of this state. 

Higher-order correction to the differential frequency 
shift, Eq.lQJ, arises due to terms quartic in the field 
strength £q. This fourth-order contribution is expressed 
in terms of a.c. hyperpolarizability 7(0;). The expres- 
sion for 7(u;) IT^ has a complicated energy denominator 
structure exhibiting both single- and two-photon reso- 
nances. While for the ground state there are no such res- 
onances, for the 6 Pq a two-photon resonance may oc- 
cur for 6s8p and 6s8p ^Pj intermediate states. Due to 
theoretical errors in calculations of the magic frequency 
we can not reliably predict if the two-photon resonances 
would occur. Since the resonance contributions may 
dominate 7(ci^), we can not provide a reliable estimate 
of the fourth-order frequency shift. The estimate may be 
carried out as soon as the magic frequency is measured 
with sufficient resolution. As a possible indication of the 
effect on the clock frequency, we notice that for Sr 
the resulting correction to the energy levels was a few 
mHz at a trapping laser intensity of 10 kW/cm^. This 
systematic uncertainty can be controlled by studying the 
dependence of the level shift on the laser intensity p] . 

The 6 "^Pq state decays due to an admixture from J = 1 
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states caused by the hyperfine interaction. In this pa- 
per we restrict our attention to the hyperfine interac- 
tion due to the nuclear magnetic moment /i. We write 



hfs 



N 



.(1) 



where tensor 



this interaction as H] 

tI^-* acts on the electronic coordinates and /LtAr is the nu- 
clear magneton. We employ the following nuclear param- 
eters: for ^^^Yb, the nuclear spin 7 = 1/2 and magnetic 
moment iJ, = 0.4919 /xat, and for ^'''^Yb, I = 5/2 and 
= — 0.6776 /LtAT. Using first-order perturbation theory, 
the HFS-induced transition rate is given by 



A 



hfs 
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where the sum S is defined as 

(6^5o||i^||7') {i\\TP\\&^PS) 



E{i)-E{Q^PS) 



(4) 



(5) 



and ujo = E{& 3Po°) - E (e^^o) . 

To estimate the rate we restricted the summation 
over intermediate states to the nearest-energy 6 and 
e^Pf states. Using the CH-MBPT method we com- 
puted HFS couplings, {Q^Pg\\r^^\\Q^P^) = -6685MHz 
and (6 ^poj jT^^) ||6 ^PC) = 4019 MHz and wc inferred the 
values of dmole matrix elements from lifetime measure- 
ments 0, [3 . The resulting HFS-induced lifetimes of 
the 6 3Po° level are 20 and 23 seconds for i^^Yb and "^Yb 
isotopes respectively. 

A coherence of atomic states may be lost due to 
scattering of laser photons (Rayleigh and Raman pro- 
cesses ^OJ). These are second order processes. In partic- 
ular, the Rayleigh (heating) rate for both 6 ^Pq and the 
ground states may be expressed in terms of a.c. polariz- 
ability 
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where II is the intensity of laser. At the magic frequency 
uj* the values of a.c. polarizability for both states are 
equal to 160 a.u. (see Fig. For a laser intensity of 
10 kW/cm^, the resulting rate is in the order of 10"'^ 
sec~^. As to the Raman rates, there are no Raman tran- 
sitions originating from the ground state. The final states 
for transitions from the 6'^Pq are the J = 1,2 sublevels 
of the same 6'^Pj fine-structure multiplet. We estimate 
this rate by approximating the relevant second-order sum 
with the dominant contribution from the 6s7s ^Si inter- 
mediate state. The resulting Raman scattering rate is 
also in the order of 10"'^ sec^^ for 10 kW/cm^ laser. 

The total magnetic moment of the Yb atom is com- 
posed of the nuclear and electronic magnetic moments. 
Disregarding shielding of externally applied magnetic 
fields by atomic electrons, the g-factor due to the nuclear 
moment is given by ^^nuc = — (l/wp)(/x//L(Ar)//, where 
nip is the proton mass. The numerical values of Sg^^c 
are -5.4 x 10"'' for "^Yh and 1.48 x 10"'* for ^^^Yb. 



The electronic magnetic moment of the 6 ^Pq state arises 
due to mixing of levels caused by the hyperfine interac- 
tion, i.e., the same mechanism that causes the 6'^Pq state 
to decay radiatively. This correction may be expressed 
as 



SgMs 
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The computed values of the Sg^fs correction are 2.9 x 
10-^ for i^^Yb and -8.1 x 10"^ for ^^^yb, which im- 
ply that mHz shifts would be produced by fiG magnetic 
fields. Fields can readily be calibrated and stabilized to 
this level using magnetic shielding. 

The hyperfine interaction also induces residual vector 
(axial) a^p (lo) and tensor a^p {uj) a.c. polarizabilities. 
For J = levels there is no tensor contribution for the 
^^^Yb isotope (/ = 1/2) and it can be shown that for the 
173 Yb (/ — 5/2) it vanishes when the HFS interaction is 
restricted to the dominant magnetic-dipole term. For a 
non-zero degree of circular polarization A, the relevant 
correction to the light shift of level 7P is 



M 

'2F 



A a^p (lu) -So 
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where for J = 0, F — I and M is the magnetic quan- 
tum number. Using third-order perturbation theory and 
a formalism of quasi-energy states |l7l | we arrived at an 
expression for a^p (oj) which contains two dipole and one 
hyperfine operator in various orderings and double sum- 
mations over intermediate states. Analyzing these ex- 
pressions, we find that the vector polarizability of the 
6'^Pq state is much larger than that for the ground state, 
as in the case of Sr For Sr, Katori et al. Q estimated 
the vector polarizability by adding HFS correction to the 
energy levels of intermediate states in Eq. Our anal- 
ysis is more complete and we find that the dominant ef- 
fect is not due to corrections to the energy levels, but it is 
rather due to perturbation of the 6'^Pq state by the HFS 
operator. The resulting values of a^apo (uj*) are —0.10 

a.u. for ^""^Yb and 0.075 a.u. for ^""^Yb. Using these 
values in the above equation, we find that holding A to 
< 10~^ with fields of 10 kW/cm^ would keep shifts in 
the clock frequency below the mHz level. This require- 
ment on optical polarization is not an extreme one, and 
in the special case of a ID optical lattice could be relaxed 
significantly by orienting the quantization axis (defined 
by the external magnetic field) perpendicular to the trap 
axis. 

In conclusion, we have analyzed the possibility of cre- 
ating a highly precise optical clock operating on the 
6 ^5*0 — » 6 "^Pq transition in odd isotopes of atomic Yb. 
According to our calculations, the natural linewidth is 
about 10 mHz, and the magic wavelength for producing 
zero Stark shift of this transition in an optical lattice trap 
is about 752 nm. We have examined possible sources of 
shifts and broadening due to both the optical trapping 
fields and any magnetic fields, and find they should not 
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perturb the clock above the 10""'^* level, except for possi- 
ble larger near-resonant terms in the hyperpolarizability. 
An accurate measurement of the magic wavelength will 
be needed to settle this last question. 
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